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FORWARD 


The research results contained in this technical report were performed under the 
NASA grant NAG-1-1684 entitled "Experimental and Numerical Structural Acoustic 
Control for Interior Noise Reduction". The report is based essentially on partial progress 
of the Ph. D. dissertation prepared by Jeffrey S. Bevan under direct guidance of Dr. Chuh 
Mei. The document presents a finite element formulation and control of sound radiated 
from cylindrical panels embedded with piezoceramic actuators. The extended MIN6 
shallow shell element is fully electrical-structural coupled. A piezoelectric modal 
actuator participation (PMAP) is defined which indicates the actuator performance to 
each of the offending modes. Genetic algorithm is also employed to validate the sensor 
and actuator locations determined by the PMAP criteria. The work was conducted at the 
Department of Aerospace Engineering, Old Dominion University. Mr. Travis L. Turner, 
Structural Acoustics Branch, NASA Langley Research Center is the technical monitor. 
The research is continued under the new grant NAG-1 -2 141. 
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1. Shallow Shell Finite Element Formulation 


The three-node shallow shell element presented is an extension of the MIN6 developed 
by Tessler 1 . A brief description of the element attributes includes 3-node shallow curved 
shell triangle with fifteen nodal degrees-of-freedom (dof). The element employs C° 
anisoparametric interpolation to account for both membrane and transverse nodal 
displacement. The membrane and transverse displacement each utilize unique 
interpolation polynomials, which offer distinctly different polynomial degrees as 
introduced by Tessler 2 . The transverse displacement employs a complete quadratic 
polynomial, while the bending rotations utilizes a linear interpolation scheme. The MIN6 
formulation consists of a Marguerre shallow shell membrane strain consistent with the 
Reissner-Mindlin theory which includes transverse shear effects. Static condensation 
eliminates the extra nodes required to facilitate the complete quadratic polynomials while 
maintaining the minimum specified number of nodal dof. By formulating a co-rotational 
highly nonlinear solution, Barut 3 et al exhibits the versatility of the MIN6 shallow shell 
element. This research further expands the MIN6 capability by including anisotropic 
piezoceramic materials in conjunction with unsymetric laminated composite resulting in a 
fully electrical-structural coupled shallow shell finite element formulation. 

2. Element Displacement Functions 

Displacement field components u x , u y , and u w , consistent with Mindlin theory, are 
described as 


u x = u{x,y,t) + zG y {x,y,t) 
ii y =v{x,y,t) + zO x {x,y,t) 
= w(x,y,t) 


( 2 . 1 ) 




where u, v, w represent the mid-surface membrane and transverse displacements; bending 
rotations of the normal about the x and y axes are given by 6 X and d v respectively. The 
element coordinate system is defined as positive x along the side between nodes 1-2 as 
shown in Figure 1. The arbitrary shallow shell shape is described by h 0 (x,y) and related 
to the r-axis as 
z=z-h 0 (x,y) 




Figure 1 MIN6 Geometric description 




Thus, the nodal displacement vectors can be written as 


M 7 ={L w *J 

KJ KJ} 

(2.2) 

M 7 =K 

W 2 K 

(2.3) 


9 ,2 0,1 0y2 0,1 J 

(2.4) 

K} 7 =h 

M, » 3 V, V, J 

(2.5) 


where each electrical dof is the coupled, total electric potential of each piezoceramic 
layer. For example, consider np layers of piezoceramic material the electric potential dof 
is given by 

W r =b, - kJ (2.6) 

Considering the inherent electrical-mechanical coupling of piezoceramics, the electrical 
potential dof includes both self-generated, or sensor voltage, and externally applied, or 
actuation voltage. The coupled piezoelectric constitutive relation is described in Sect. 4. 


The displacement field throughout the element is determined by interpolating the 
nodal displacement as 


'*(*. >’.') = L # . JH } + L H.S M 

(2.7) 

= Lft * ftj IM+U. ft L, M, Aft 

bft,M=Lft * ft o o oJ{e) 

(2.8) 

s,(*,.y,<)=L/ft>}=[p 0 0 ft ft ftojjs) 

(2.9) 

Ii(x,y,t)= LW,J{w.}=Lft ft ft 0 0 °J 

(2.10) 

v(.v,y,/)=Lw.JK}=L» 0 0 ft ft ftj 

(2.11) 
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Natural coordinates commonly used to describe triangles refer simply to area 


ratios. These area or natural coordinates are related to the geometric 

coordinates by utilizing the following transformation relations 
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X 2 

X 2 

^2 * 
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( 2 . 12 ) 
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where (x, , _y ( ) designate the i' h nodal coordinate, and the triangular area A is given by 
A = ((x 2 - x, )(y 3 -_V|)-(x 3 - x, )(y 2 - y \ )). The interpolation functions are defined as 


follows 


L 2 = '-(b,N 5 -b,N t ) 

i, = '-(!>, 2N, -ft, AT,) M, = jM, -a,N.) 

M, = jM, -a,N,) M,=j(a,N, -a,N t ) 

W 4 = (2.13) 

a, = x 3 , a, = x /3 a 3 - x,, 

b, = y i2 b 2 = y l} b, = y 2i 

x ii = x , ~ x j y u 

Integration of polynomials expressed in area coordinates can be accomplished with the 
following formula 


\£ k £,4 m dA=2A 


kUlml 

(2 + k + l + rnf. 


(2.14) 
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Geometric matrices describing the strain interpolation functions are given by the 
following definitions. 
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3. Strain-Displacement Relations 

Tessler 1 first developed the shallow shell element utilized by combining both 
Reissner-Mindlin and Marguerre theories. The Marguerre theory enforces the Kirchhoff 
assumption thereby neglecting shear deformation in contrast to Reissner-Mindlin theory. 
Substituting the normal bending rotations for the overall slope in the Reissner-Mindlin 
theory yields the following Marguerre inplane strain compatible to Reissner-Mindlin 




r 

-sf 
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> — < 

K’y 0 , I 
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(3.1) 


The following curvature and shear strain completes the strain relations 




(3.3) 
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4. Constitutive Relations 

The k >h layer of the laminate specifies either structural or piezoceramic properties and 
is characterized by the following coupled constitutive relations 


H=bl(M-£»{4) o.i) 

M=[alW (4.2) 


A* - {4 _£?]* (f } 


K 


{A}=-[ fi dW — 


0 



(4.3) 


(4.4) 


Eq. (4.3) defines D ik the tf h layer electric displacement density resulting from two 
components, namely strain and the total electric field Ejk, which are coupled through the 
piezoelectric d coefficients. The piezoceramic layer is polarized in the 5-direction and 
anisotropic in the 1-2 and 1-3 material directions. Furthermore, each piezoceramic layer 
may have an arbitrary orientation angle. The lamina reduced stiffness components are 
determined from the principal material properties as 


, — VA, — — 

1 _ V \2 V 2\ ' " 1 — V 12 v 2l (4.5) 

Qc, 6 =G n Q u =G 23 0 55 = G l2 

The ability to accurately model piezoceramic anisotropy supports current research 
trends in advanced transducer development. If anisotropic piezoceramic material is used 
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the following transformation relates the principal coordinate charge coefficients to the 
global coordinates 


'd x 
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cos~a 
sin 1 a 
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c/ 32 > 
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(4.6) 


5. Force and Moment Resultant 

Analysis of laminated composite plates maintains distinct lamina stress, therefore 
utilizing stress resultants is imperative. The stress resultants, or force and moment per 
unit length are defined as 

({A (UK (5.1) 

/2 k 



dz 


Utilizing Eq.(5.1) it is useful to define the stress resultants as follows 


(5.2) 


H- 

[M Mil 

W 

Lfo [r>lJl 



(5.3) 


where the extension, extension-bending, and bending stiffness matrices are defined as 


M (**+■ -**) 

k = I k 

Z 1 k 

[D)4±[Q] (£,-$) 

J *=l k 


(5.4) 

(5.5) 

(5.6) 

(5.7) 


Considering the k! h piezoceramic layer and the coupled constitutive relations the 


piezoceramic force and moment vectors are given by 
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(5.8) 


({A'4-H))=j|[eL{< £ .. ('•?)£ 

6. Equations of Motion 

Finite element equations of motion for the laminated composite panel with fully 
coupled electrical-structural properties are derived utilizing the generalized Hamilton's 
principle 4 to obtain 

f 2 8(J-U + W, -W m +W)dt = 0 (6.1) 

■M 

where T and U are the kinetic energy and strain energy of the system, W t is the electrical 
energy, W m is the magnetic energy, and W is the work done due to external forces and 
applied electric field. The magnetic energy is negligible for piezoceramic materials if no 
external magnetic fields are located near the specimen. The kinetic energy of the shallow 
shell element is defined as 

r=j)p({w} r W + { 1 i}’{ 1 i} + {v} r {v}^ (6.2) 

V- “ 

where w, u , and v are the transverse and membrane velocity components and p is the 
mass per unit volume, and — is the volume of the element. The potential and electrical 
energy are defined as 

(6.3) 

V- “ 

W' - \ X -{E} r {D}dP (6.4) 

H - 

and the work done on the element by external sources is defined as 
W = J {w}‘ {F h }dV + J{w} r {F f )dS + {w} r {F c } - \VpJS (6.5) 

* Sx S i 
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where {F a } is the body force vector, {F} is the surface traction vector. {F t .} is the 
concentrated loading vector, S ] is the surface area of the applied traction, S 2 is the surface 
area of the piezoelectric material, V is the voltage applied to the piezoceramic layer, and 
p cs is the total electrical charge due to self-generated piezoelectricity and applied 
actuation voltage. In Hamilton's principle, all variations must vanish at the time t = r, 
and t = t 2 . The Hamilton's variational statement may be written in the most general form 
as 

jMw' {M + faY {«} + {&}' W) 

- {&}>} + {«}'{£>} + {Sw? {F t )W 

+ j(&v)' 1^, 'f K - \syp„dS * W IF 1=0 (6.6) 

s, Si 

Evaluation of Eq. (6.6) leads to the development of the finite element matrices and the 
elemental equations of motion. Employing the stress resultants, the variational potential 
and electrical energy may be described as 

S(U-W t )= \[{Se"} r {N} + {5 k}'' {M} + a {8y} r {R}- {5 Ef {D}yA (6.7) 

A 

where the shear correction factor for the MIN6 element is defined as 



The finite element equations can be determined by completing the variational work 
statement in terms of the nodal values. Furthermore, it is convenient to specify the 
shallow shell slope matrix in the following form 


k]= 


K„, o 

0 K, y 

h<) ’V ’ X 


(6.9) 
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where h 0 =h 0 (x,y) and is the constant shallow shell geometry. The stress resultants are 
written in nodal components, Eq. (5.5) as follows 


M=Mk}+[olW-K} 

=[ B ][c.]{ w .}-[ fl lK]{«}+[o][c,]{tf}-fA/ # } 


Variation of the electrical energy requires attention due to z dependence 




a L 


dz 


dA 


Completing the integration with respect to z yields 


( 6 . 10 ) 

( 6 . 11 ) 

( 6 . 12 ) 


(6.13) 


ftp _ i 

A _ ^-1 ^ 

-(^.){<[e], E »K +(<5£«)C£ m h t ]]dA 


(6.14) 


Before continuing, we can further expand the definitions of the piezoelectric force and 
moment vectors. The force vector may be expressed as 


M=I 

* = 1 

- KLkU 
= -[[§], M,*, - [elK*. 

=-kMKl 




QlMKK [fid!"'.) 


Similarly, the piezoelectric moment vector can be expressed as 


(6.15) 
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(6.16) 


where 



Eq.(6. 1 4) can be recast in matrix form using the above definitions as 

= J{^} r ([^vfK}+[/’»,]'W-[)r]K]‘'{£ 1 ) + [^][B,]‘ l {£ J })^ ( 6 . 18 ) 

where 

rwIelM - o 

H= : Uu\o\U\< i (6.19) 

0 - 

Thus, the expansion of the variational potential and electrical energy can be achieved 
by substituting the aforementioned relations as follows 

S(u-W)= {jV} + {&}' {M} + a {Syf {«} - {SE) r {D}Uf 

= j{({^.l r [c.]' -W K]')(M[C-„]{vn.}-MK]{^}+[ J B][C-»]{^}-{)V,}) 

A 

+ a({5vv # } r [c rf ] , +{^} r [c^] r )([4][c„]{w,} + [4][c^]{(?}) 

+ {^)'[ s d'[[ ;, vnc.]{w.} + [p v ] r [e , ,]{<?}+[/=; ,nc.]{s} 

+ {[ e »]-[z]){ w i}_ dA 

( 6 . 20 ) 

Expanding Eq.(6.20) term by term yields expressions leading to the element stiffness 
matrices 


1 1 




= j{l^dWM[c.lK} 

A 

(6.21) 


(6.22) 


(6.23) 


(6.24) 


(6.25) 


(6.26) 


(6.27) 

+ W1«J[^1«,K} 

(6.28) 

+ {«} r [c,J'[sIcJ w .} 

(6.29) 

-WlcJMf.M 

(6.30) 

*{se}‘\c b ][D\c h le) 

(6.31) 

+wi c h \p u \B t \« t \ 

(6.32) 


(6.33) 

+ a {6»j[cJ[A,\cJ{0} 

(6.34) 

+a{seY[cJ[A,\c»lw t } 

(6.35) 

+a{sey{cj[A,\cje } 

(6.36) 

- I**',!' ]' [ p n T \c. Ifw, } 

(6.37) 

-W'WMfcM 

(6.38) 

-Kfkl'hJlc.M 

(6.39) 

+ !<>“, }' \b, ]' fe l-hBk I}'- 1 

(6.40) 


Completing the generalized Hamilton's principle considering nodal degrees-of- 
freedom yields inertia, external mechanical loading and piezoceramic actuation quantities 
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^ r (6.41) 

-P^™,,,} 1 {{ H „}[ H „ J { K, } - { J { w„ } )} dA - j{Sw,} {p a }dA 

s r 

Finite element stiffness matrices result from examination of the potential and electrical 
energies of the variational work statement. Each stiffness matrix, including fully coupled 
electrical-structural and geometrical shallow shell stiffness’ is defined from the indicated 


term of Eqs.(6.2 1 )-(6.4 1 ) 

{56}' [k 0 ]{#} where [k e ] = J[C„ ]' [D\C h J iA (6.42) 

A 

where [*J= j[C»]>I ’C.ytA (6.43) 

A 

{ Bw „, \ k „,6 M where [k m0 ] = J[C m ]' [b\P„ \iA (6.44) 

A 

1' l k „, ]K, } where [k m ] = j[C„, ] ? [a\C„, \iA (6.45) 

A 

\ k „,4 JK } Where [Kt ] = /I Cm \P N \ B 0 V 4 (6-46) 

A 

}' [ k <„, ] K, } where [k m } = J[5, ]' [P v ]' [C m \lA (6.47) 

geometric stiffness due to shallow shell geometry 

{fa * }' l k n,o L {°] where [ k me L = JI C ,„ r U\ d „ \ iA (6- 48 ) 

A 

{S0\‘ [ k a„ L K, } where [k a „ ], = J[<9„ ]' [A\C m )lA (6.49) 

A 

I se) r [k,\\e) where [k,l . Jtf,n^K]-t!»J[BlC,I-[C 4 n«Ie.j)« (6.50) 

A 
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WM.K) where KA = 

(6.51) 

K where K1 = \M[ p siMA 

A 

(6.52) 

stiffness due to shear effects 

a i s °y KM where [*J= WcJ'UslcM 

A 

(6.53) 

a { 5d Y[ k M»} where kj= HcjuMW 

A 

(6.54) 

Y k 0 ]W where K, ] = Jfo ]' [A 5 \C yd \iA 

A 

(6.55) 

a { Sw * Y K } where K ] = j[c^ ] 7 [A s \lA 

A 

(6.56) 

coupled piezoelectric-structural stiffness 

W [ k e* IK I where [k (H ] = J[c* f [P m \b* \lA 

A 

(6.57) 

I<H f Ko IM where [k„ ] = J[^ J 7 [P M J 7 [C h \lA 

A 

(6.58) 

Kv )' K IK 1 where [fcj= {K J (K" 3 ]- [/]K 

(6.59) 

mass matrices and load vectors from Eq. (6.41) 

k]= 

A 

(6.60) 

["> K ] = \[HJ p[H«]dA 

A 

(6.61) 

["•mh \{HJ p[H.]dA 

A 

(6.62) 

[>",]= \\ H J p[ h A m 

(6.63) 
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(6.64) 


K]=/([^.] + [«.])V([W.]+[ff.])<M 

A 

{. P b ( 0 } = p(x,y,t)dA (6.65) 


{Pv(‘)} = p(x,y,t)dA 

A 


( 6 . 66 ) 


[pA t )} = ~\PcM (6.67) 

A 

Resulting in the following finite element equation of motion 


K] [ m M>] 0 

[' m OH ] K] 0 

0 


/ 

''0 

0 

0 

0 

0 

9 

> + 

0 

M+m,, 

[b„] + [b„l 

[b] + [b]„ 

0 0 0 

K] 

< 

0 

[Ke]-[K, 6 \ 

[*-] 

[b*] 

0 0 0 

0 

, V V 

V 

0 

[*# ]-[**]. 

[bn] 

M 


k ] [*, 1 0 0 

\ J L J 

\ 

r 

w h 


A (O' 

.*r 

i i 

i 1 

1 1 

O 

o 


e 

> = < 

Pa(0 

0 0 0 0 
0 0 0 0 

\) 

W ,n 

.V 


0 

p* (Oj 


( 6 . 68 ) 


Through standard finite element assembly lead to 


M 0 
0 0 





Separation of Eq.(6.69) produces two coupled equations of motions 
[Mp} + ([£. ] - [K, t I K, ]-' [K„ ]M = {/>. } - [k,, \K , ]" {/> } 


(6.69) 


(6.70) 

(6.71) 
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Eq.(6.70) describes the system with respect to the primary variable {W}, the structural 
nodal dof, and the applied loads, indicated by the {/(,.}, Subsequently, 

Eq.(6.71) defines the electrical dof as a function both the structural dof and the external 
actuation voltage. Thus, finite element formulation supports simultaneous sensing and 
actuation. However, if the piezoceramic patches are restricted to actuation or sensing 
only, the fully coupled formulation is preserved. 

7. Modal Formulation 

The number of nodal degrees of freedom can be excessively large and thus 
impede the controller design and computational efficiency. However, the modal 
coordinate transformation can be utilized by considering the overall structural response 
can be very well represented by a reduced number of normal modes. Thus, expressing 
Eqs. (6.70) and (6.71) in the modal domain yields 

Wr M + [C r fe}+ [K, ] = {F „ } - \F r ^ } (7.1) 

K}=k]"h}-k]"KI^fe} (7.2) 

where; [m r ] = [y/J MM , [K, ] = [yrj ([*C„ ] - [k w$ \k< ]'' [k^ \ys] , { F nv } = [///] 7 [P w ) , 

and } = W \ K * ] ' {^» }■ The damping matrix [C r ] will be determined 
experimentally from the corresponding damping ratios of the individual modes measured. 
The modal matrix, or eigenvectors, [>//] represents a reduced set of the first ten normal 
modes. Hence, the large number of finite element equations simply reduces to ten 
equations in the modal domain. 
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8. Development of Actuator Location Optimization Process 

The piezoelectric modal actuator participation (PMAP) [A/]”'[^] 7 J -1 {/> (*)} 

indicates the contribution of the actuators to each of the modes. Exploiting the coupled 
electrical-structural formulation and the modal domain properties the PMAP is evaluated 
and ranked for desired acoustical modal contribution. Hence, the preferred actuator 
location is determined. Results determined using a rectangular 24-dof plate finite 
element with three piezoceramic patches are shown in Table 1. The patch numbering 
refers to Figure 2, where patch 2 is indicated by a 2, and patches 8 and 9 are shown by 8 
and 9 respectively. Note that patch 9 consist of nine finite elements and patches 2 and 8 
each have one element. Furthermore, each acoustic radiating mode is indicated along 
with the corresponding PMAP. The PMAP values indicated in Table 1 are normalized 
indicating the maximum actuator performance for mode (1,3) with patch 9. The 
piezoceramic stress values </j;=rf?.r=0.171e-9 was used and is representative of PZT-5A 5 . 


Patch 2 

Patch 8 

Patch 9 

Mode 


0.2151 

0.2151 

-6.8134 

(U) 


0.0901 

0.0901 

- 1.0000 

(1,3) 


0.0562 

0.0562 

0.2868 

(3,1) 


-0.0592 

-0.0592 

0.0072 

(5,1) 


-0.0747 

-0.0747 

-0.2661 

(3,3) 


-0.1326 

0.1326 

- 0.0000 

(4,3) 



Table 1 PAMP for 10”xl4' , x0.040” rectangular plate 
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nap2 = 


11111112 1 
111111111 
1119 9 9 111 

1119 9 9 111 

1 119 9 9 1 1 1 

111111111 
181111111 ] 

Figure 2 Vector representing typical piezoceramic sensor-actuator location map of three 
patches. Key: (1) no piezoceramic present, (2) a single patch, (8) a single patch, and (9) 

patch comprised of nine finite elements. 

9. COUPLED ACOUSTICS USING RADIATION FILTERS 
The structural velocity distribution may be constructed by superposition of a set of 
independent acoustic modal velocities^ 1 . The acoustic velocity modes correspond to the 
acoustic pressure distributions of the radiated sound of the structure. Thus, a radiation 
filter may be constructed by considering a radiation operator which represents 
independent discrete acoustic structural radiators. The radiation operator is frequency 
dependent and may be subjected to a singular value decomposition, which yields 
interesting and useful results. Consider the following relationship of the total radiated 
power of a vibrating structure 

{ P } = K/'*>F [*(/®)M/<»)} (9- 1 ) 

where {v(yo>)} is a nxl vector of nodal velocities and [/?(y<y)] is a nxn radiation matrix. 
The radiation matrix is proportional to the radiation resistance with diagonal and off 
diagonal elements corresponding to self and mutual radiation resistance respectively. 
Thus, a vibrating structure represented by n discrete independent acoustic point sources 
or radiators represented by the following radiation matrix 
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M= 


to p,,A‘ 


4/TC 


1 


sin 


fan) 


sin fai2) . sinfa,,) 

sin(Ar 23 ) 


At-,, 

sin fa„i) 


At* 12 

I 


kr, 


I n 


kr„ 


kr ' 23 

sin fa;»-.) 


(9.2) 


where co is the circular frequency, p Q is the density of air, A is the elemental area of the 
corresponding radiator, c the speed of sound in air, k the wave number (o/c), and r pq is 
the distance between the p and velocity location. Singular value decomposition of 
the radiation matrix yields dominant acoustic radiation modes for each frequency. Thus 

[*(/<»)) = [■»]'' (9.3) 

where [«] = [{«}, {«}, {«}„] is a nxn matrix whose column represent normalized 


acoustic radiation modes at co p and [l] is a nxn matrix of singular values The 

magnitude of the q' h singular value indicates the relative importance of the corresponding 
acoustic radiation mode. The singular values, as a function of frequency, represent the 
dominant coupled acoustic modes, which provides guidance for developing a noise 
control strategy. The associated normalized singular vectors obtained from the singular 
value decomposition represent the contributing acoustic radiation mode. Figure 4 shows 
the dominate acoustic radiation modes determined by evaluating the radiation matrix 
singular values. Furthermore, Figure 5 presents the dominant corresponding acoustical 
SVD radiation modes. 


19 








;ue i 


Mode 2 



Mode 3 




Mode 5 



Figure 6 Acoustic Radiation SVD Mode Shapes 




10. Control Strategy 


The objective of this research emphasizes determining the preferred placement of 
piezoceramic sensors and actuators locations to minimizing the structure-borne radiated 
noise of a cylindrical shell, to this end J~C 2 and LQG 7 control strategies will be employed. 

The J-fy and LQG has been well researched and applied to laboratory experiments. 
Although it may not represent the best practical control approach much data has been 
published and will provide guidance. The state space approach is achievable by 
considering the following fully coupled equations of motion 

[A ip} + i[K w )- [K.,}[K t ]){W} 

Using a modal transformation results in the following set of uncoupled modal equations 


f 

q r +2<Z r Q) r q r +co;q r =-i- 


( 10 . 2 ) 


Where the modal mass and stiffness are obtained from 


and 

/, =vi f ({^(0}-fcJL > l' , !n(^) (io.4) 

The following state space equations for LQG can be utilized 

j4)}=.[4K/)} + [Bj U (,)} + [B.M0} (10.5) 
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(>’(')} = [C WO} + {v(/)} , {x} = {q qf 

where the following state space transformation was applied to the finite element coupled 
equations of motion 



[ 0 ] 

WrV[K r ] 



r i T [o] 1 r [o] 

[8 ‘ = [\m, ]-' WX [k„ \k, ]-' |J = [[M r ]- [^ f _ 

{Ch\K,Y[K^\ [ 0 ]] 


Evaluation of the system corresponding to the previous piezoceramic patches was 
conducted by observing the open loop impulse response as shown in Figure 7. Each 
patch is shown in Figure 2 and designated by 2,8, and 9. Again, the data represents a 
structure modeled using rectangular plate elements modified to include the couple 
electrical-structural formulation. 
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